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REMOTE PROFILING OF LAKE ICE USING AN S-BAND SHORT-PULSE 
PvADAR ABOARD AN ALL- TERRAIN VEHICLE 
by Dale W. Cooper, Robert A. Mueller, and Ronald J. Schertler 
NASA Lewis Research Center, Cleveland, Ohio 44135 

ABSTRACT 

An airborne short-pulse radar system to measure ice thickness 
was designed and operated during the 1973-74 and 1974-75 ice seasons. 

This system supported a joint effort among NASA, the U. S. Coast Guard 
and NOAA to develop an all-weather Great Lakes Ic,r ■ \formation System 
which aids in extending the winter navigation soa?' This paper describes 
experimental studies into the accuracy and limitatiohi of this short-pulse 
radar system. A low power version of this system was built and operated 
from an all-terrain vehicle on the Straits of Mackinac during March 1975. 
The vehicle allowed rapid surveying of large areas and eliminated the 

I 

ambiguity in location between the radar system and the "ground truth" 
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ice auger team which is unavoidable for the airborne versions of the system. 
These in situ measurements also permitted an assessment of the effects of 
snow cover, surface melt water, pressure ridging, and ice type upon th? 
accuracy of the system. Over 25 sites were explored which had ice thick- 
ness in the range of 29 to 60 cm. The maximum radar overestimate was 
9.8 percent, while the maximum underestimate was 6. 6 percent. The 
average error of the 25 measurements was 0. 1 percent. 

Introduction 

During the past three years NASA, in conjunction with the U. S. 
Coast Guard, National Weather Service, National Oceanographic and 
Atmospheric Administration, and U.S. Army, has been developing an 
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all-weather Great Lakes ice information system to aid in extending the 
winter navigation s^son. An extended season has the potential of saving 
millions of dollars in coal and ore shipping costs, since cargo must now be 
shipped by more costly rail or truck routes or stored until spring thaws. 

The entir>i operational ice information system is scheduled to be turned 
over to the U.S. Coast Guard by the end of the 1975-76 ice season. 

The Great Lakes ice information system primarily uses a side- 
looking airborne radar (SLAR) aboard a U.S. Coast Guard C-130 aircraft. 
This X-band radar is an improved version of the system flown on Army 
OV-1 aircraft in Viet Nam. The SLAR provides an all-weather .•'erial view 
of the ice cover which is than transmitted to SMS (synchronous meteoro- 
logical t ’Ute). From the satellite the information is retransmitted to 
a NOAA .. on at Wallops Island, and then via telephone lines to the U. S. 
Coast Guard's Great Lakes Ice Center at Cleveland for interpretation. 
Detailed ice maps are then constructed and sent by radio facsimile to any 
ship on the Great Lakes 'with the appropriate recording equipment. This 
information is updated daily to reflect ice shifts due to wind. 

The SLAR system has proved to be A^ery effective in determining 
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ice location patterns and movement on the Great Lakes in virtually all types 
of weather. The SLAR is very sensitive to surface roughness and discon- 
tinuties and readily gives the location of pressure ridges. Because the 
surface pattern is most often a relic of the early history of the ice, the 
SLAR imagery cannot be interpreted directly to give ice thickness. To 
supplement the SLAR data ice auger teams have been used in the past, but 
these ground truth measurements are laborious, time consuming, expensive, 
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weather dependent, dangerous to personnel and cannot be done on a large 
enough scale to map an entire lake in a reasonable amount of time. It is 
for this reason that a short pulse (one nanosecond) S-band radar system 
was developed to profile the thickness of ice remotely. 

The remote ice measuring system was fi it designed and 
tested for feasibility during the winter of 1972-73 and flown successfully 
aboard a U.S. Coast Guard H-53 helicopter at altitudes up to 100 meters 
[Vickers et al. , 1973]. Other radar methods were considered [Page et aL , 
1973: lizuka et al., 1971], but appeared to have possible operational and 
resolution problems. 

For the next ice season (1973-74) the nanosecond radar system 
was redesigned to be used as an ice profiler aboard a C-47 aircraft 
[Cooper et al. , 1974]. This system "'roved operational at altitudes up to 
2300 meters and ground speeds of 75 m/sec. The radar was able to detect 
ice thickness from 10 cm to 92 cm, the upper limit being the thickest ice 
found on the Great Lakes in the past two ice seasons. 

As thickness data from the short-pulse radar system was 
utilized in an operational program in support of the SLAR during ice sea- 
sons 1973-74 and 1974-75, questions arose as to the radar's accuracy and 
operational limitations. During some initial helicopter check out flights 
a small amount of calibration data was taken. Some more calibration of 
the short-pulse radar system on the C-47 aircraft was done in conjunction 
with an ice auger team on Brevoort Lake in the upper oeninsula of Michigan 
west of the Straits of Mackinac during the ice season of 1973-74. The lake 
was of relatively uniform thickness and easy to locate by aircraft, i The 
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radar data checked the auger team to within 2 cm, Because K was not 
possible to exactly locate the ice auger team on a flight line, some uncer- 
tainty in the calibration remained. 

Operationally, results from the C-47 flights had shown that the 
radar was unable to detect less than 10 cm of ice and that measurements 
were precluded wherever surface melt water exceeded about 1 mm in 
thickness, because of lack of penetration. Snow cover on the ice was never 
found to be a problem unless the snow surface was slushy. The S-band radar 
system worked well in any type of weather except rain which, of course, 
has the same effect as melt water. 

I 

Calibration of the system could not be performed readily in the 
laboratory because of the lai^e amount of ice and water required to approxi- 
mate a planar surface. In the laboratory it also would be difficult to make 
any supporting structure nonreflective to microwaves. The dielectric con- 
stant of ice made from pure water has been measured in many laboratory 
studies [Birks, 1961; Gumming, 1952; Evans, 19G5], Von Hippie [1954] 
found the dielectric constant of ice to be 3,2 at S-band, decreasing to 3. 17 
at X-band. Accurate measurements of ice thickness depend on a good 
knowledge of the dielectric constant of lake ice. 

To help facilitfct > better radar measurements, ten lake ice sam- 
ples were taken from the Straits of Mackinac and Whitefish Bay during the 
winter of 1973-73 and analyzed for dielectric constant and loss tangent at 
Stanford Research Institute [Vickers, 1975]. The diameter of the samples 
was 7. 6 cm With thickness ranging from 30 to 78 cm. Various ice types 
were in the samples, i. e. , clear, milky with small air inclusions, and 
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moi;dy clear with large layered inclusions. Ice machining to permit the 
insertion of the sample in the test wave guide was performed in an environ- 
mental chamber at -10° C. The relative dielectric constant for lake ice 
from 1 to 18GHz was found to be 3. 17 for clear ice with no inclusions, 3.08 
for milky ice with small (less than 5x10"^ cm diameter) air inclusions and 
2. 99 for clear ice with large (greater than 0. 6 cm) air inclusions. Whereas 
these results are in the expected range of values, the laboratory method 
requires ice storage for long periods of time prior to examining as well as 
machining which could possibly alter the measured results. 

Because locating an ice auger team directly on a flight line for 
accurate system calibration was considered nearly impossible, it was 
decided to build a low power radar which could be mounted on an all-terrain 
yehiclr . With this configuration the radar could be easily checked without 
location ambiguity. A low r.f. power radar was required to insure person- 
nel safety. A portable gasoline generator supplied the auxiliary 115v, 60Hz; 
electrical power. The all-terrain vehicle was borrowed from the Coastal 
Zone Labs of the University of Michigan at Traverse City, Michigan. The 
actual research using the all-terrain vehicle on the Straits of Mackinur was 
performed in the second week of March. 

!< THEORY AND SYSTEM CONFIGURATION 

The short-puise radar system works on the basic radar principle 
of time delay. From an aircraft the radar pulse may be considered a plane 
wave which is partially reflected upon incidence at the surface of the air-ice 
or snow-ice interface. Part of this wave continues on through the ice at a 
slower velocity due to the increased relative dielectric constant. At the 
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ice-water interface total reflection takes place and the wave continues at 
its slow rate back to the air-ice or snow-ice interface^ The receiver 
measures the time delay between the partially reflected incident wave 
and the wave which continued through the ice to the water interface and 
back before preceedir^ to the aircrafts To know the velocity of the wave 
in the ice, the relative dielectric constant must be known. In addition to 
ice thickness, a measure of the time that the partially reflected incident 
wave takes to return to the aircraft gives aircraft altitude. 

The time delay t in nanoseconds can be directly related to the 
ice thickness, x in cm, by the following expression if the relative dielec- 
tric constant of ice, is known: 

X 14.99 t 

Generally a relative dielectric constant of 3. 1 [Vickers, 1975] is assumed 
for lake ice which yields: 

X 8. 51 t ■ 

js 

The short-pulse radar system has linade Ice thickness measure- 
ments in support of the SLAR imagery during the ice seasons of 1973-74 
and 1974-75. As previously mentioned, aircraft location is of prime impor- 
tance because radar makes measurements at the nadir the aircraft. NASA's 
C-47 aircraft is equipped with an inertial navigation system (INS) which is 
accurate to ± 1. 85 Km af^;er one hour. While the INS is adequate for 
making ice charts, good calibration cannot be performed. 

Both C and S-band versions have been employed on the C-47 air- 
craft, but the S-band proved operationally superior because of better system 
components. The aircraft radar system block diagram is shown in Figure 1 


7 


The S-band system used , 4ther random noise or continuous wave modula- 
tion at 2. 86 GHz. Random noise modulation was used to avoid the possi- 
bility of coherent interference between the transmitted pulse and other 
interfering signals. In actual operation this problem did not occur as 
theorized. 

The heart of the system is the nanosecond pulse generator 
which when mixed with the S-band oscillator allows only a few cycles of 
r.f. power to be transmitted. For this purpose, a double mixer system 
is enqployed to decrease feed thru from each double balanced mixer. The 
length of coaxial tr'«jismission line to the second mixer is cut to the proper 
length so both pulses,, the output of the first mixer and the output of the 
pulse generator, arrive simultaneously. The traveling wave tube ampli- 
fier gives over 20 watts of peak power. The solid state pre-amplifier 
maintains a proper input level. The entire system' bandwidth must be 
greater than 1 GHz, as dictated by the nanosecond pulse. l 

For the purpose of narrowing the receiver antenna pattern, 4 
ridged horns feed a combiner and then a low noise (3. 8 db noise figure) 
solid state amplifier, all located in close proximity. After detection a 
1 GHz sampling oscilloscope was used for the final display. The radar 
was initially triggered by the clock which operated from 40 to 250 kHz. 

The oscilliscope was triggered from a delay unit which started the scope 
at the precise time that the return pulses were received, A manual 
adjustment was used on the delay unit which required constait manipula- 
tion by the operator as the aircraft altitude changed to kei^'p the data dis- 
played on the oscilloscope face. Recording of data was done with an 
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oscilloscope camera. Typical aircraft ice return displays liave been 
reported [Cooper et al., 1974 Vickers et al. , 1973], A new design is 
now being formulated which determines the ice thickness electronically 
and corrects for any aircraft altitude changes. In addition this new system 
will profile the ice surface, permitting the detection of ice ridges. 

ALL-TERRAIN VEHICLE CONFIGURATION 

When applying radar principles to an all-terrain vehicle, the 
nearness of the ice involved new considerations. It was necessary to de- 
termine (hat the ice surface was in the far, not near, field of the transmit- 
ting antenna. Near field calculations were made on the S-band pyramidal 
ridged horn [Harrington, 1961], At least 20 cm of separation were required 
between the transmitting horn and the ice surface. 

Personnel safety was of primary in^ortance because of long 

exposure times, so the radio frequency power density in any personnel areas 

2 

was maintained below 1 mW/cm . This was a full decade below usual U.S. 
safe standards. Fortunately this still allowed sufficient power to make 
thickness measurements. The 115V, 60 Hz power was supplied from the 
auxiliary gasoline generator. 

The radar for an all-terrain vehicle is simpler than the airborne 
system because of the lower radio frequency power levels involved. The 
all-terrain vehicle S-band short pulse radar block diagram is shown in 
Figure 2. This system differs from Figure 1 in that there was no random 
noise source, transmitter traveling wave tube amplifier and delay unit. 

The solid state amplifier was the final transmitter with only 10 dBm out- 
i^jiut which insured personnel safety. Triggering of the sampling oscillo- 
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scope was directly from the clock without a delay unit. 

The system was composed physically of four subsystems, i, e. , 
the clock, the sampling oscilloscope, thfj iransmitter and receiving radar 
module and the two ridged pyramidal horns. The transmitting and receiv- 
ing modulo with power supplies weighed 2.7 Kg. 

The all-terrain vehicJe in operation is shown is Figure 3. This 
vehicle allowed 3 persons to be carried at speeds up to 48 km/hr. Due to 
the low bottom of the vehicle, deep snow and high ridges were impassable, ji 

Polyethylene bags were placed over the horn antenna,s to keep mud, snow, 
and water out of them in transit. With clean bags over the antennas the 
system was con^letely operational on the ice. The antennas were extended 
1.5 m from the rear deck of the ail-terrain vehicle and both were adjusted 
in the traverse plane while pointing downward to have maximum output. The 
transmitting horn’s source of radiation wa,s taken to be 1.30 m above the ice 
surface, while the receiving horns point of collection was taken to be 1. 14 m 
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above the ice surface. The distance between these points measured parallel 
to the ice was 1. 25 m. 

In reducing the data from the all-terrain vehicle, it was neces- 
sary to take into account the geometry of the situation. Since the receiving 
and transhiitting horns were in a nonsymmetrical configuration the geomet- 
rical theory of diffraction [Kline et al. , 1965] was employed in data reduc- 
tion. Computer iterations using ray tracing techniques at various thick- 
nesses accounted for refraction at the ice surface. The wave propagation 
was entirely normal in the aircraft model. 


The wave bending as it enters the ice may be obtained from Snell’s 
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law for electromagnetic waves [Harrington, 1961]. The angle to the normal 
in air at the ice surface, 9, may be related to the angle to the normal in ice, 
(p, by the relative dielectric constant of ice, 

For this ail- terrain vehicle radar system, the computerized ray tracing 
technique yielded the following approximate equation for ice thickness over 
10 cm: 

X « 8,77 t if = 8.1 

r ;i 

A series solution can also be obtained [Cooper et al. , to be published 1976] 
that yields the same result. Where ice is not identifiable, a dielectric 
constant of ice, of 3. 1 was assumed. 

On the C '47 aircraft a refracted wave, directly from the transmit- 
ter to receiver, was an inconsequential problem because essentially it 
occurred microseconds prior to the two return pulses separated only by 
nanoseconds. It could readily be removed from the receiver return. 

On the other hand, in the all-terrain vehicle configuration due to 
the short distances, this refracted pulse exists in all the data and there- 
fore must be disregarded. The first pulse is the refracted pulse, the sec- 
ond pulse is a surface reflected wave, while the third pulse is reflected 
from the water as shown in Figure 4. 

Data were originally to be taken with a polaroid oscilloscope 
camera. The severe cold caused the shutter to stick and visual data could 
only be taken with a 35 mm camera. A typical ice return is shown in 
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Figure 5. This photograph ‘«t^as taken on the Straits of Mackinac near 
St. Helena where the auger team measured 39 cm. The feed-thru pulse 
occurs first followed by the pulse from the ice surface and finally by the 
pulse from the ice-water interface. The measured radar delay is 4.2 nan- 
oseconds or 36.8 cm of ice if the relative dielectric constant is 3.1. This 
ice was covered by 8 cm of snow which is not detectable. The icu appeared 
milky with some similar air inclusions. 

TEST RESULTS 

The short pulse radar Was calibrated by comparing actual auger 
and radar measurements. This calibration was done on the Straits of 
\^ckinac from March 12-14, 19'‘/5. 

Over twenty-five sites were examined and measurements, were 
compared. Considerable effort was made to find varied sites, i.e. , different 
ice types and thicknesses, Ocassionally deep snow made the all-terrain 
vehiclf; inoperable. Essentially all the encountered ice was snow covered 
as shown in Figure 3. Up to 25 cm of snow was found at one site but it 
did not affect any of the radar measurements except when the snow had a 
slushy crust which was only found at one location. Snow was usually re- 
moved around the site to help identify the ice type. Large blocks were not 
removed as had been in previous ice studies at the Straits of Mackinac. 
Usually a dielectric constant of 3. 1 was assumed. 

Previously the airborne radar sometimes became inoperable in 
warm temperatures. Surface melting was theorized as the problem. Exper- 
iments were performed which verified the fact that melt water prevents pene- 
tration of the radar wave into the ice. Tests showed that de-ionized water 
was as effective as lake water in preventing penetration. 
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The all-terrain vehicle was very limited in crossing ridges due 
to its low underside. Because there was fear of becoming immobile in deep 
snow, hourly check-ins by marine radio were made to the local U. S. Coast 
Guard. 

Near the Mackinac Bridge, a pressure ridge was studied where 
two slabs of ice had refrozen, one over the other. The radar quickly re- 
vealed an ice step from 29 cm to 50 cm. Auger measurements supported 
this finding, A high rafted area 82 cm thick gave inaccurate radar results 
of about 35 cm but investigation revealed that the ice had a slushy interior. 
Another rafted area measured 110 cm with the auger, but the radar gave a 
reading of only 93 cm. Large air pockets were present indicating that the 
relative dielectrir ;'\?;'Stant, was assumed too high at 3, 1. 

For uniform ice with small air inclusions the radar and auger 
meaBuremenis were in close agreement assuming t j. = 3 . 1 in over 25 differ- 
ent sites as shown in Figure 6. The line represents exact agreement be- 
tween auger and radar measurements. These sites had ice thickness in the 
range from 20 to 60 cm. 

DISCUSSION OF ERROR 

Some of the data scatter in Figure 6 was directly attributable 
to the ice auger team. The ice auger measurements could be up to ±1 cm 
in error due to tape reading, a local convexity or concavity at the bore hole 
and top and/or bottom surface roughness. This estimate has taken into 
account considerable past ground truth experience which has shown the bot- 
tom surface of lake ice to be relatively smooth and flat. 
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The radar measurements shown in Figure 6 were normalized to 
the auger measurements for comparison. The maximum radar overestimate 
was 9. 8 percent, while the maximum underestimate was C.6 percent. The 
average error of the 25 measurements was 0.1 percent. 

A major unknown in the radar mesurements is the true dielectric 
constant of the ice. In fact, all measurements could be considered exact if 
the dielectric constant of the ice varied from 2. 8 to 3.8. As can be seen by 
the magnitude o! these numbers, they are out of range of what ice is thought 
to be and other factors are contributing to the errors. 

In the present configuration the pulse time delay can only be 
accurately read to about 0.2 nanosecond which corresponds to a 1.75 cm 

.. il_ 

thicknesb error which is about 5%pf a 35 cm measurement. 

Another concern is the location of the actual, electrical interface. 
It could be above the actual ice-water interface or the ice and v/ater could 
have an air pocket between them. 

The aircx'aft system has not been able to measure ice below about 
10 cm because it is not possible to separate the two return pulses. This 
limit has not been verified with the all-terrain vehicle radar system because 
no ice was found at the proper thickness. A measurement of 92 cm of lake 
ice has been detected from an aircraft, which is the thickest lake ice we have 
found in the past two ice seasons. 

CONCLUSION 

The short pulse radar system has been demonstrated to be an 
accurate remote lake ice measuring device. Snow cover or adverse weather 
generally did not interfere with making measurements. Surface melting 
and rain however preclude measurements. 
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Further studies are necessary to establish the variability of 
the relative dielectric constant for lake ice. The best meam^ to do this 

'i 

is in situ, possibly with the short pulse radar and an auger t^iam. For 

next ice season a radar system is being built for the Coast Guard's 

I 

HV-16 helicopter stationed at Traverse City, Michigan. T!iis could allow 
many different types and ages of ice to be studied throughout the ice season. 

NASA’s OV-1 and Coast Guard's C-130 aircraft will have C-band 
versions of this radar system for next winter. At the higher frequency, 
smaller antennas will be used, but surface roughness may be more of a 
problem. A half nanosecond pulser will be used to increase the resolution 

of the instrument. 

1 '.. 

Electronic processing of the data will be done and manual setting 
of the oscilloscope trigger delay will no Icjiger be required. Pressure 
ridges will be studied b^J noting any rapid changes in the first return pulse. 
Knowing the pressure ridge height relative to the surrounding ice field 
should yield ridge thickness, assuming that the ridge is in isostatic balance, 
i. e. , its weight is supported by its bouyancy, not by stresses in the sur- 
rounding ice field. 
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